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Abstract 

We propose a computational procedure for creating a stable equilibrium triple junction (TJ) with con- 
trolled grain misorientations. We apply this procedure to construct a TJ between a S5(210) grain boundary 
(GB) and two general high-angle GBs in copper, and calculate the diffusion coefficients along the TJ and the 
GBs using molecular dynamics with an embedded-atom potential. The TJ diffusion is only a factor of two 
faster than diffusion in the S5 GB but significantly faster than diffusion in the general GBs. Both the GBs 
and the TJ studied here show a premelting behavior near the bulk melting point, where their diffusivities 
converge to the diffusivity of bulk liquid. Although our results are consistent with the common assumption 
that TJ diffusion is generally faster than GB diffusion, the difference between the two diffusivities does not 
appear to be large enough to ensure a significant contribution of TJs to diffusional creep in polycrystals at 
high temperatures. 
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I. INTRODUCTION 



The triple junction (TJ) is a line where three grain boundaries (GBs) meet together. It is com- 
monly believed that TJs greatly influence properties of nanocrystalline materials and can play a 
significant role during plastic deformation and diffusional creep!^ 

The TJs are much less studied crystalline defects than GBs or dislocations are. Experimental in- 
vestigations of TJs are complicated by many factors, including their high mobility."^ Experimental 
measurements indicate that zinc diffusion along TJs in Al is three orders of magnitude faster than 
along GBs."^ TJs can also be the sites of preferred solute segregation.^ Some theoretical models 
introduce the assumption of fast TJ diffusion to explain the anomalous diffusion and deforma- 
tion behavior of nanocrystalline materials. However, no reliable experimental information on 
self-diffusion along TJs is currently available to verify this assumption. 

Due to the recent progress in the development of interatomic potentials capturing basic prop- 
erties of materials on the quantitative level,'^! atomistic simulations offer a means of predicting 
self-diffusion coefficients that would otherwise be very difficult to measure. For example, the 
GB self-diffusion coefficients in copper predicted by simulations were found to be in very good 
agreement with experiment.!^ To the best of our knowledge, TJ diffusion has never been studied 
by atomistic simulations. 

As in experiments, the TJ motion at finite temperatures poses a major challenge to diffusion 
simulations. In a few studies published so far, the TJ motion was prevented by creating high- 
symmetry configuration^ or restricting the simulations to the K temperature. Applying these 
restrictions, the atomic structure of the TJs at K was carefully examined and their excess energy 
was estimated/^^^ however, no diffusion coefficients were computed. TJ motion by capillary 
forces at high temperatures has been studied by two-dimensional Lennard- Jones simulations.^^ 

In this work we apply three-dimensional molecular dynamics (MD) simulations with an 
embedded-atom potential to compute self-diffusion coefficients along a representative TJ and to 
compare them with self-diffusion coefficient in the adjoining GBs. Since self-diffusion is the only 
type of diffusion studied here, we will be referring to it for brevity as simply "diffusion". In Sec- 
tion [ll] we show how a TJ with controlled grain misorientations can be created in the computer and 
how its motion can be constrained by boundary conditions. We also introduce our methodology 



for diffusion calculations. Our results are presented in Section HI followed by their discussion 
and conclusions in Section ITVl 
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II. METHODOLOGY 



We chose copper as the model material, with atomic interactions described by the embedded- 
atom potential developed in Ref. [141 The same potential was used in the recent study of GB 
diffusiorP'. The bulk melting point of Cu predicted by this potential is Tm = 1327 K (the experi- 
mental value is 1356 K). 

Our procedure for constructing an equilibrium TJ is illustrated in Figure [T] First, a sym- 
metrical tilt S5(210) boundary was created with the tilt axis parallel to the common [100] di- 
rection of both grains. The grain dimensions are 85 x 83 x 66 A, where the third number 
is the dimension normal to the page. The GB was terminated at a free surface perpendicular to 
the tilt axis. Another 116 x 66 x 66 A grain with the orientation shown in Fig. [l] was placed on 
the surface, forming a non-equilibrium TJ. The model was equilibrated by a 2 ns long MD run at 
1 100 K with a periodic boundary condition parallel to the TJ and fixed boundary conditions at the 
bottom and the sides of the block (the top side is exposed to vacuum). The fixed layers, shown in 
Fig.[l]in gray, were 10 A thick. Further details of the MD methodology will be discussed later. 

During the MD anneal, the TJ quickly (within the first 0. 1 ns) moved down, causing migration 
of the two initially horizontal GBs until an equilibrium configuration was reached. The planes of 
the moving GBs changed in this process, but the grain misorientation remained the same. The new 
GB planes were symmetrical with respect to the S5 GB by construction, but they did not align 
with any particular low-index crystal planes in the grains. These new boundaries will be referred 
to as general GBs (GGBs). 

Once equilibrated, the TJ and the GBs did not move except for small thermal fluctuations. 
Note that the S5 GB is pinned by the fixed region at the bottom while the GGBs are pinned by the 
surface grooves formed in the comers between the third grain and the two initial grains. 

Finally, a cylindrical region with the axis parallel to the TJ was cut out of the tricrystal 
(Fig. [T|c)). The atoms within a 10 A thick outer shell of the cylinder were made fixed, whereas 
all atoms inside the shell were allowed to move during the subsequent MD simulations. The fixed 
shell served to pin the GBs and prevent their migration. Note that the cross-section of the shell 
is not perfectly circular: it contains the previously formed surface grooves and a portion of the 
bottom fixed layer from the original block. The boundary condition parallel to the TJ remained 
periodic. This cylindrical block contained 81,137 atoms and had the length of 66 A. It was used 
as the initial configurations for the diffusion calculations. 
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The MD simulations were performed in the NVT (canonical) ensemble using a Nose-Hoover 
thermostat. The integration of the equations of motion employed the Verlet algorithm with the in- 
tegration time step of 2 fs. The simulations covered the temperature range from 700 K (the lowest 
temperature of reliable diffusion calculations accessible with our computational resources) to 1315 
K (12 K below T^). Prior to an MD run, the block was uniformly expanded/contracted according 
to the thermal expansion factor at the intended simulation temperature, using the thermal expan- 
sion factors computed in Ref . TT, The model was then additionally equilibrated by a 1 ns long MD 
run to approach the equilibrium point defect concentrations at the chosen temperature, followed by 
a 10 ns long production run for diffusion calculations. As a test of equilibration, the total energy 
of the simulation block was averaged over several time intervals covering the production stage, 
and it was checked that these average values followed a normal distribution. Similar tests were 
made for the pressure inside the block. An additional proof of sufficient relaxation comes from 
the observations that (i) during the subsequent diffusion runs (see below), the TJ displacements 
were random and very small (< 5 A), and (ii) the mean- squared atomic displacements accurately 
followed the Einstein relation. Snapshots containing the atomic coordinates were saved every 0.02 
ns and later post-processed to compute the diffusion coefficients. 

The diffusion coefficients D were calculated from the linear fits of {X"^) versus At plots using 
the Einstein relation D = (X^)/2At . Here, At is the diffusion time and (X^) is the mean-squared 
displacement of atoms over the time At in the X direction parallel to the TJ. Since the atoms are 
labeled, (X^) was computed by comparing the initial and final coordinates of all atoms initially 
contained within a selected probe region and averaging their squared displacements in X over 
the time interval between t and t + At. This obtained value of (X^) was additionally averaged 
over all choices of the reference time t by increasing it from zero to tmax — At, where tmax is the 
duration of the entire MD run (10 ns). The maximum diffusion time At was chosen to be 1 to 2 
ns, depending on the temperature. 

The probe regions selected for GB diffusion were parallelepipeds aligned with the GB planes 
(Fig.|2|, which had the dimensions LxW x 2d, where L = 66 A is the cylinder length. The width 
W = 35 A parallel to the GB was chosen so that to include as large a GB area as possible but 
avoid the influence of the TJ and the fixed regions. Since our goal was to evaluate the diffusivity in 
the GB core area, the half-thickness of the probe region was d = 5 A, which is close to the typical 
diffusion width of GBs.ti^l 

For TJ diffusion, the probe region was a cylinder with the radius of d placed at the diffusion 
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center of the TJ (Fig. |2]). Special care was taken to keep track of the position of the center. To 
achieve this, a testing MD run was made prior to the production run at each temperature and a 
group of atoms with the largest magnitude of displacements in X was identified. Such "fastest" 
atoms were invariably found in the TJ region, which was already a qualitative proof of its enhanced 
diffusivity. The average coordinates of the "fastest" atoms were used to select the diffusion center 
of the TJ. Furthermore, the position of the center was recomputed for each snapshot during the 
production simulations and it was found to change only by small amounts that did not exceed d. 
This confirmed the hight stability of our TJ geometry. However, the position of the center had to 
be slightly re-adjusted for each temperature since changes in the contact angles of the GBs could 
produce TJ displacements. 

Diffusion coefficients in liquid copper were computed separately from mean-squared displace- 
ments of atoms at several temperatures using a large periodic supercell at zero pressure. 

III. RESULTS 

The TJ and GB structures were examined after "quenching" from 700 K to K. Although 
the structures contained a number of quenched-in point defects, the were very coherent and had 
thicknesses on the order of 2-3 interatomic distances. These structures remained quite coherent 
also at finite, but not too high, temperatures. From the dihedral angle of the GBs, the K energy 
of the GGB is estimated to be approximately 0.9 of the S5 GB energy. This number is very 
reasonable, given that the S5(210) GB has the highest energy among all boundaries studied in 
Refs. ll6ll"7l and is roughly 10% above the energies of typical high-angle GBs. 

Examination of the MD snapshots reveals that the GBs and the TJ become increasingly dis- 
ordered with temperature, as illustrated in Fig. |3] (The visualization employs the coordination- 
number mode of the AtomEye program,^^ plotting only atoms whose coordination number is 
smaller than a certain number.) Near the bulk melting point the GBs become increasingly wider, 
especially in the TJ region, and develop a liquid-like structure. At 1315 K the TJ essentially turns 
into a liquid pipe. 

Typical (X^) versus time plots are displayed in Fig. Q. Their linearity demonstrates that the 
Einstein relation is followed accurately, confirming that we have properly sampled the diffusion 
regime. The slopes of the curves indicate that the TJ diffusion is faster than diffusion in the S5 
GB, and that both diffusivities increase with temperature. 
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All diffusion coefficients obtained in this work are summarized in the Arrhenius diagram in 
Fig. |5] Observe the encouraging agreement with experimental data for self-diffusion in the bulk 
liquid^^ and in high-purity polycrystalline Cup^j^ lending confidence to our simulation method- 
ology. Diffusion along the TJ is seen to be faster than diffusion in both GBs at all temperatures 
tested. Also, the S5(210) GB shows a higher diffusivity than the GOB at all temperatures. The 
latter observation is consistent with the previous work"', in which diffusion in the S5(210) GB 
was found to be faster than in the E5(310) and E17(530) GBs. Note the GB and TJ diffusivities 
converge at high temperatures and tend the diffusivity of bulk liquid, which is in agreement with 
the premelting process. To put these data in perspective, the self-diffusion coefficient in the lattice 
computed with the same EAM potential at the bulk melting point is as low as 1.1 x 10^^'^ m^/s, 
which is beyond the scale of Fig. |5| 

Despite the significant scatter of the data points, the Arrhenius plots in Fig. [5] can be linearized 
except near the melting point. A fit according to the Arrhenius equation D = DqCxp {—Q/kT) 
(k is the Boltzmann factor) in the temperature interval 700-1100 K gives the activation energies 
Q and pre-exponential factors Dq summarized in Table |lj Although these Arrhenius parameters 
slightly depend on the choice of the temperature interval, their ordering remains the same. The 
activation energy of TJ diffusion is always the smallest and close to that of the S5 GB. 

IV. DISCUSSION AND CONCLUSIONS 

We have proposed a computational method to create a single equilibrium TJ with controlled 
crystallographic orientations of the grains. The geometric stability of the TJ is ensured by the 
fixed boundary conditions which pin the GBs at three points. An alternative procedure could be 
to prepare three separate wedge-shaped grains with angles that sum to 180°and join them together 
like a cake. Besides being much more tedious, this approach does not offer an easy way to elimi- 
nate the mismatch stresses that may arise due to possible changes in the translational states of the 
grains during the relaxation. In addition, the free volumes of the GBs and of the TJ would have to 
be accommodated by elastic stresses. In contrast, in our method the TJ forms naturally, as a result 
of GB migration at a high temperature (Fig.[TJb)). In this process, the equilibrium atomic density 
in the GBs and in the TJ can be reached by atomic diffusion to/from the free surface. Although 
in this paper we have implemented a particular set of grain misorientations, shown in Fig.[T] other 
grain misorientations could have been created as well. A precise control the GB planes is more 
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difficult, although some control is possible through judicious choices of the dimensions of the 
initial grains. 

The MD simulations indicate that both the GBs and the TJ become increasingly disordered 
at high temperatures and exhibit a premelting behavior near the bulk melting point (Fig. |3]). As 
the temperature approaches T^, the premelting first starts at the TJ, then extends to the S5 GB, 
and finally occurs at the GGBs. This sequence of premelting processes suggests that the TJ has 
a higher excess free energy than the GBs studied here. The question of whether TJs carry an 
excess free energy relative to the adjoining GBs ("excess of excesses") has long been debated in 
the literaturcfi^ 

Diffusion simulations have been conducted over a wide temperature range and the diffusion 
coefficients in the TJ and the GBs have been extracted from the mean-squared atomic displace- 
ments. It should be mentioned that due to the fixed boundary conditions and conserved number 
of atoms in the model, the point-defect equilibrium at the diffusion temperature could only be 
achieved by defect generation and/or absorption at the TJ and the GBs. At high temperatures, the 
TJ and the boundaries were premelted or at least highly disordered, making them effective sinks 
and sources of point defects. Even at the lowest temperature studied here (700 K), the TJ region 
was still fairly disordered and was presumably able to generate and absorb point defects during the 
pre-diffusion anneal. This was evidenced by the linearity of the mean-square displacements versus 
time plots at this temperature, as well as by the independence of the diffusion results on the ther- 
mal history. For example, an additional low-temperature anneal of the TJ prior to carving out the 
cylinder (i.e., when the GBs were still connected to the surface, Fig.[TJb)) did not change the diffu- 
sion coefficients within the statistical error. It is important to recognize, however, low-temperature 
MD simulations in closed systems containing of well-ordered GBs or other extended defects can 
be associated with non-equilibrium point-defect concentrations and can give inaccurate diffusion 
coefficients. 

Our calculations confirm that diffusion along the TJ is faster than diffusion along the GBs at 
all temperatures (Fig.|5]). Our results also suggest that, for diffusion along "typical" GBs, such as 
the GGB tested here, can be orders of magnitude slower than TJ diffusion at low temperatures. 
This conclusion is consistent with recent experimental measurements.^ However, diffusion in the 
highly energetic S5(210) GB is only a factor of two slower than in the TJ. Although this small 
gap tends to slightly increase at low temperatures. Fig. |5] and the Arrhenius parameters (Table |I]) 
clearly indicate that high-diffusivity GBs such as S5 can be more efficient for atomic transport in 
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polycrystals than TJs. The factor of two higher diffusivity of TJs is hardly sufficient for making 
a significant contribution considering their much smaller cross-sectional area in comparison with 
GBs. 

As the temperature approaches the melting point, the TJ and GB diffusivities converge to the 
bulk liquid diffusivity, reflecting the premelting process. A very similar behavior was recently ob- 
served for GB diffusion-fiSl The rapid increase of the diffusivity near Tm produces a strong upward 
curvature of the Arrhenius plots and the concept of activation energy loses its physical signifi- 
cance. From the convergence of the diffusivities at high temperatures, it appears unlikely that TJs 
can make a sensible contribution to diffusional creep at high temperatures, given their relatively 
small cross-sectional area even in nanocrystalline materials. With most GBs in nanocrystals being 
of general type, the contribution to the creep rates often attributed to TJs might be overestimated. 
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Table I: Arrhenius parameters Q (activation energy) and Dq (pre-exponential factor) of TJ and GB diffusion 
in Cu obtained in this work. The Arrhenius parameters of lattice self-diffusion computed with the same 
interatomic potential are given for comparison. 



Defect 


Q(eV) 


Dq (m2/s) 


Triple junction 


0.47 ±0.02 


(4.3 ± 1.3) X 10"^ 


S5GB 


0.48 ±0.04 


(2.5 ± 1.4) X 10"^ 


General GB 


0.71 ±0.02 


(1.0 ±0.2) X 10""^ 


Perfect lattice 


i.96r 


3.1 X 10^6^ 



"Sum of the vacancy formation and migration energies (Table 2 in Ref. [T6l) 

''Computed from the lattice parameter and the vacancy formation entropy and jump-attempt frequency 
(Table 2 in Ref. [Hi 




Figure 1: Schematic showing the procedure for the creation of a TJ. (a) Three grains are joined together by 
putting a third grain on top of a bicrystal. (b) An equilibrium TJ forms during an MD run at a high tempera- 
ture, (c) A cylindrical region containing the TJ is cut out of the tricrystal and additionally equilibrated. The 
shaded areas represent fixed atoms. 
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Figure 2: The plot indicates the probe regions chosen for diffusion calculations for the TJ and the GBs. 
Only atoms with coordination number below a certain bound are shown. 
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700 K 1100 K 1315 K 

Figure 3: Change of the structure of the GBs and the TJ with temperature. The visuaUzation method is 
explained in the text. The bulk melting temperature with this interatomic potential is 1327 K. 
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Figure 5: Arrhenius diagram of diffusion coefficients D calculated for tiie TJ (A), the S5 GB (□), tlie GGB 
(O) and liquid copper (V). The vertical dashed line indicates the melting point of Cu with this embedded- 
atom potential (1327 K). Experimental data for self-diffusion in liquid coppei^^ (+) and in high-purity 
copper polycrystals^" ( x ) are shown for comparison. 
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